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Spatial structure of an individual Mn acceptor in GaAs 



O 

o 

(N 



C/2 



A.M. Yakunin,! A.Yu. Silov,i P.M. Koenraad,i J.H. Woltcr,i 
W. Van Roy,2 J. De Boeck,^ J.-M. Tang,^ and M.E. Flatte^ 

'cobra Inter- University, Eindhoven University of Technology, 
P.O.Box 513, NL- 5600MB Eindhoven, The Netherlands 
^IMEC, Kapeldreef 75, B-3001 Leuven, Belgium 
^Optical Science and Technology Center and Department of Physics 
and Astronomy, University of Iowa, Iowa City, lA 52242, USA 
(Dated: February 2, 2008) 

The wave function of a hole bound to an individual Mn acceptor in GaAs is spatially mapped by 
scanning tunneling microscopy at room temperature and an anisotropic, cross-like shape is observed. 
The spatial structure is compared with that from an envelope-function, effective mass model, and 
from a tight-binding model. This demonstrates that anisotropy arising from the cubic symmetry 
of the GaAs crystal produces the cross-like shape for the hole wave-function. Thus the coupling 
between Mn dopants in GaMnAs mediated by such holes will be highly anisotropic. 

PACS numbers: 71.55.Eq, 73.20.-r, 75.50.Pp 
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Despite intense study of deep acceptors in III-V semi- 
conductors such as Mnca, little information has been ob- 
tained on their electronic properties at the atomic scale. 
Yet the spatial shape of the Mn acceptor state will in- 
fluence hole-mediated Mn-Mn coupling and thus all of 
the magnetic properties of hole-mediated ferromagnetic 
semiconductors 1] such as Gai-a-Mn^jAs. Evidence for 
anisotropic spatial structure even in shallow acceptors 
such as Zn, Cd 0, 1^ and C Q in various III-V semi- 
conductors suggest that anisotropic hole states may be 
common. In addition to controlling inter-dopant proper- 
ties such as the Mn-Mn coupling , the wave- function 
shape would affect single-dopant properties such as the 
g-factor and optical transition oscillator strengths. 

This study presents an experimental and theoretical 
description of the spatial symmetry of the Mn acceptor 
wave-function in GaAs, and we suggest our results im- 
ply similar behavior for other acceptors and other hosts. 
We first present our measurements of the spatial map- 
ping of the anisotropic wave function of a hole localized 
at a Mn acceptor. To achieve this, we have used the 
STM tip not only to image the Mn acceptor but also to 
manipulate its charge state /A~ at room temperature 
as described in 7]. Within an envelope function effec- 
tive mass model (EFM) the anisotropy in the acceptor 
wave-function can be traced to differing amplitudes of 
envelope functions with the same total angular momen- 
tum (L > 0) but different angular momentum projec- 
tions along a fixed axis. We introduce into the EFM a 
single parameter r/ that describes the breaking, by the cu- 
bic crystal, of spherical symmetry for the acceptor level 
envelope functions. As 77 has a negligible effect on the 
binding energy compared to the central cell correction, 
common variational approaches cannot be used to eval- 
uate rj. However, comparison with calculations based 
on a tight-binding model (TBM) for the Mn acceptor 
structure jQj permits us to clearly identify the physical 
origin of the anisotropic shape in these models, and to 
justify the value of rj used in our EFM calculations to 



describe the experimental shape. The TBM calculations 
also demonstrate that although the spin-orbit interaction 
does influence the acceptor wave- function, the qualitative 
anisotropic shape of the acceptor state occurs in crystals 
without spin-orbit interaction. Thus acceptor levels in 
crystals such as GaN should have a similar shape. 

The measurements were performed on several samples 
using chemically etched tungsten tips. The samples con- 
sisted of a 1200 nm thick layer of GaAs doped with Mn 
at 3 X 10^* cm"'^ grown by MBE on an intrinsic (001) 
GaAs substrate. A growth temperature of 580 °C was 
chosen to prevent the appearance of structural defects 
such as As antisites, which would complicate the spatial 
mapping by shifting the position of the Fermi level of the 
sample. The concentration of the Mn dopants was low 
enough to neglect Mn-Mn interactions and the formation 
of an impurity band. The samples we used were insulat- 
ing below 77 K. The experiments were performed in a 
room temperature UHV-STM (F < 2 x 10~" torr) on 
an in situ cleavage induced (110) surface. 

A major advantage to our approach is that the occu- 
pation of the acceptor state can be influenced by band- 
bending from the voltage applied between the STM tip 
and the sample (see the Fig. ^ inset). We studied the 
voltage dependent appearance of the Mn acceptor in the 
STM constant-current image. In the ionized configura- 
tion at high negative bias Mn appeared as an isotropic 
round elevation which is a consequence of the influence 
of the A~' ion Coulomb fleld on the valence band states 
[Fig. |2Ia)]. This agrees with a recent study of the indi- 
vidual Mn in GaAs in the ionized conflguration |^]. We 
found that at a positive bias Mn is electrically neutral, as 
can be seen from the absence of the electronic contrast at 
high positive voltage {U > 1.5 V) when the conduction 
band empty states dominate tunneling (see the positive 
branch of the I{V) curve in Fig.^). At low positive volt- 
age where the tip Fermi level is below the conduction 
band edge Mn appeared as a highly anisotropic cross like 
feature [Figs.[2Ib) and|3|^a)]. The anisotropy is even more 
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FIG. 1: I(V) curves acquired on the clean GaAs surface (solid 
line) and in the middle of the cross- like feature(dashed line). 
The simulated position of the flat-band potential Ufb is indi- 
cated by the arrow. Insert displays energy band diagram for 
the positive sample-bias. 



evident in a reciprocal-space image [Fig.Efb)]. The meth- 
ods of calculating the theoretical images [Figs.[21[d),inic- 
f)] will be described below. 

The cross like feature manifested itself in the local tun- 
neling I{V) spectroscopy at low voltages when the GaAs 
bands do not contribute to the tunneling. It appeared 
as an empty states or filled states current channel in the 
band gap of GaAs depending on applied positive or neg- 
ative bias, respectively. Thus the mapping of the Mn 
acceptor state in the filled (empty) states mode was re- 
alized by electron (hole) injection into the A'^ /A^ state. 
In the tunneling I(V) spectroscopy the manganese 
channel appeared presumably above the flat band po- 
tential Ufb and was available for tunneling in the wide 
range of voltages above Ufb- Our estimated value of 
Ufb is about -1-0.6 V. The observed ionization energy, 
which was determined from the shift of the I{V) spec- 
trum at negative bias, corresponds to the Mn acceptor 
binding energy Ea = 0.1 eV. 

The concentration of the dopants we observed with 
STM corresponds to the intentional 3 x 10^^ cm~^ dop- 
ing level. All of the dopants could be found either in the 
ionized A~ or the neutral A'^ charge state depending on 
either negative or positive sample bias respectively. In 
the experiment we identified Mn located in at least 6 dif- 
ferent layers under the surface. In order to determine the 
actual position of the Mn dopants we analyzed the inten- 
sity of the electronic contrast of the Mn related features. 
Based on the symmetry of the cross like feature superim- 
posed on the surface lattice we distinguished whether the 
dopant is located in an even or odd sub-surface layer. We 
found that at any depth the shape of the hole on {110} 
surface had nearly C21, symmetry around the surface nor- 
mal, and was elongated along the [001] direction relative 
to the [110] direction. The cross-like features induced by 
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FIG. 2: (a) 8x8 nm^ STM image of an ionized Mn acquired 
at -0.7 V; (b) 5.6x5 nm^ STM image of neutral Mn acquired 
at -1-0.6 V. Big and small round features correspond to As 
and Ga related surface states, respectively. Presumably, Mn 
is located in 3'^'' sub-surface atomic layer; (c) A model of 
the (llO)-surface (top view) representing the Mn on Ga site 
located in an odd sub-surface atomic layer (counting surface 
layer as zero); (d) 5.6x5 nm^ simulated image (logarithm of 
local density of states) of the Mn located in the 5"* sub-surface 
atomic layer (TBM). 



deeper dopants were more elongated in the [001] direc- 
tion. The feature was weakly asymmetric with respect 
to [110] direction (lowering the symmetry to the single 
(110) mirror plane) which may come from the symmetry 
properties of the (110) surface. 

A four-band envelope-function effective Luttinger- 
Kohn Hamiltonian provided one framework (EFM) in 
which to analyze the spatial structure of the bulk- 
like neutral acceptor complex formed by a valence hole 
loosely bound to a negatively charged Mn^+3d^ core 
(Mn^+3c?^ -I- hole complex). The ground state of this 
acceptor in a zincblende semiconductor can be approx- 
imated as four-fold degenerate with a total momentum 
of the valence hole F = 3/2 and has the symmetry of 
the top of the valence band Fg 0. We neglected pos- 
sible effects caused by the presence of the (110) surface 
and quantum spin effects from the exchange interaction 
between the MnSd^ core and the hole. We also ignored 
the excited states, as the energy separation between the 
ground state i?a(l<5'3/2) = 113 meV and the first excited 
state £^0(25*3/2) = 25 meV exceeds room temperature 

According to Ref. P the acceptor wave-function in 
zincblende semiconductor is represented as a four- 
component column written in the basis of Bloch func- 
tions of the valence band top Fg. The form of the wave- 
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where c^/c^^ = ?7, < r/ < 1, (2) 

yL.m{d,^p) are spherical harmonics, i?L(r) are radial 
parts of the envelope functions, and c^^ are constants 
described below. 

As pointed out by Kohn and Schechter 0, in 
this model i?i(r)Yj, (^) are eigenfunctions of parity. 
Thus the wave function [Eq. (Q] contains only even s- 
like {L = 0) and d-like {L — 2) components. In Eq.^the 
angular part of the d-component includes components 
which transform according to the and represen- 
tations of the tetrahedral point group. Their correspond- 
ing coefRcients are the constants and c^^ , respectively, 
whose ratio is denoted here as 77 [Eq. [2] . The value of 77 
depends on the valence band parameters of the particu- 
lar material and could be evaluated variationally from the 
binding energy 0, ■ However the introduction of 77 
only weakly affects the binding energy and is negligible 
compared to a central cell correction. To clarify the na- 
ture of the Mn acceptor ground state here we used 7/ as a 
free parameter, and thus varied the spherical anisotropy 
from T] — 1 , corresponding to a completely spherically 
symmetric solution, till 77 — > , where the E component 
vanished. When 7/ = 1, the function [Eq. transforms 
to that one obtained by Baldereschi and Lipari in the 
spherical approximation '13', . 

Although the EFM describes the symmetry of the Mn 
acceptor, it fails to predict the ground-state binding en- 
ergy Eo without a central cell correction. Deep acceptors 
such as Mn in GaAs are usually described by the zero- 
range potential model 0, 0| . In this model the spatial 
extension of the wave function depends on a parameter 
a, determined by Eq. 

" = h ^ , 1 h /3= — , (3) 

n \ (i + l J nih 

where the parameter a characterizes the spatial local- 
ization of the hole trapped at the acceptor and /3 is the 
ratio of light to heavy hole masses. For GaAs we used 
a = 1.197 nm~^ and /? = 0.132. This model allows one 
to obtain analytical expressions for Ro (r) and R2 (r) . 
The admixture of the d-like component to the ground 
state calculated by Baldereschi and Lipari for GaAs 
is 30% and corresponds to what we have found using the 
zero-range potential model. As can be seen from Ref. 
the d-like component dominates at distances R > 1/a. 
The admixture of the d-like component and thus the spa- 
tial structure of the confined hole are strongly dependent 
on (3. The admixture of higher orbital momentum will 
vanish in the limit f3 ^ 1. In the case /? = the envelope 
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FIG. 3: (a) Constant-current image of the neutral Mn ac- 
quired at 4-0.9 V. Mn is presumably located in the 2""^ sub- 
surface atomic layer; (b) Fourier spectrum of the image (a); 
(c) simulated image of the Mn acceptor ground state 153/2 
(EFM); (d) Fourier spectrum of the image (c); (e) simulated 
image of the Mn acceptor DOS (TBM); (f) Fourier spectrum 
of the image (e). Images (c) and (e) were simulated assuming 
Mn to be located in the 4*'' sub-surface atomic layer. 



functions Rq and R2 can be derived in an implicit form 
[1^ . In Ref. authors emphasize the different behavior 
of i?o and R2 at large distances. The Rq drops exponen- 
tially whereas R2 has power law decay at infinity. 

In Fig. |3Jc) we present the cross-section of the bulk- 
like charge density of the total wave function in the limit 
case 77 = cut through an imaginary (110) plane at the 
distance of 0.8 nm from the Mn position. The plot is 
presented with a logarithmic scale because of the inverse 
exponential dependence of the tunnel current on the tip- 
sample distance during the actual STM experiment. Fig- 
ure Ol^d) presents the calculated Fourier transform. Note 
the presence of the satellite harmonics that arise from the 
steep fall off of the wave function in the [001] direction. 
The model reproduces the key symmetry elements in the 
Fourier spectra: the elongation of the spectra along [110] 
direction, the satellites in [001] direction and the cross 
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like shape. However the cross-hke shape and the intensity 
of the sateUites are more pronounced in the experimental 
Fourier spectrum. Based on the analyses we performed 
we conclude that only the harmonic l2,i gives rise to the 
intensity of the satellites in the Fourier spectrum. The 
cross feature in the Fourier spectrum is also influenced by 
the envelope function R2 and is more pronounced when 
i?2 has inverse power dependence on distance rather then 
exponential decay. 

The tight-binding model (TBM) we use is based 
on the deep level model of Vogl and Baranowski [l8l| . 
The dangling sp'^-bonds from the nearest-neighbor As 
hybridize with the Mn d-states of Fis character. The 
antibonding combination of these becomes the Mn ac- 
ceptor state. Coupling to the d-states of character is 
weak, and hence neglected. The hybridization strength 
is fully determined by the acceptor level energy. This 
model, if further approximated within the EFM, predicts 
77 = 0, similar to what we found by fitting the EFM to 
the experimental measurements. 

The calculations of the local density of states (LDOS) 
based on the TBM are shown in logarithmic scale in 
Fig. [2Id) and Figs. |2Ie,f). The results show symmetry 
under reflection in the (110) plane, and asymmetry un- 
der reflection in the (001) plane. For Mn dopants sev- 
eral layers down from the surface, as in Fig. |2Ib), the 
shape of the acceptor state does not depend that sen- 
sitively on the spin orientation of the Mn-core Sd-spin. 
For these dopants the shape does not depend on the spin- 
orbit interaction; we confirmed this by obtaining a similar 
cross- like acceptor structure using a tight-binding Hamil- 
tonian without spin-orbit interaction, and with empirical 
parameters designed for optimal agreement with the bulk 
band structure of GaAs p^] . The situation differs greatly 
for Mn near the surface, where the axis of extension of 
the acceptor state rotates with the spin orientation of 



the Mn core spin. At these temperatures we can expect 
the Mn core spin to point in a random direction, and for 
the STM measurements to average over the possible spin 
orientations. Thus for Figs.EJd) and|Si;e,f) the LDOS is 
averaged over the Mn spin orientation. 

Although the symmetry is well retained for all simu- 
lated positions of Mn under the surface, in both models 
the best fits were achieved when the apparent depth of 
the Mn is assumed to be 2 atomic layers larger. The 
reason for this could be the vacuum barrier, which will 
tend to shift the wave function of the Mn acceptor state 
deeper into the crystal than one would predict for the 
sliced bulk crystal calculation. The lateral size of the 
wave function measured is also somewhat larger than cal- 
culated for Fig. 131 this might come from a reduction of 
the acceptor binding energy very near the surface. 

In conclusion, we have experimentally demonstrated 
that the Mn acceptor ground state has highly anisotropic 
spatial structure. This spatial anisotropy is due to a 
significant presence of d-wave envelope functions in the 
acceptor ground state. We have demonstrated that the 
observed symmetry can be explained within simple tight- 
binding model whose only free parameter is the acceptor 
level energy. We also found that this spatial structure 
can be described well by a simple four-band envelope- 
function model of cubic symmetry whose key parameter, 
77, can be fit to the observed spatial structure, and is sim- 
ilar to that expected from a tight-binding model. These 
results have broad implications for all acceptor-acceptor 
interactions in zincblende semiconductors, and especially 
for hole-mediated ferromagnetic semiconductors. 
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